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Abstract 

There are deep similarities between non-linear QFT 
studied in high-energy and non-equilibrium physics in con- 
densed matter. Ideas such as the Schwinger mechanism and 
the VoUcov state are deeply related to non-linear transport 
and photovoltaic Hall effect in condensed matter. Here, we 
give a review on these relations. 
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Figure 1 : Several phenomena in condensed matter physics 
in strong electric fields plotted in the {E, ri)-space. 

In strong field physics, researchers are interested in the 
change of the "quantum vacuum" due to strong external 
fields. A typical example is the decay of the QED vac- 
uum in strong electric fields due to the Schwinger mecha- 
nism |[T|. When a strong enough electric field is applied 
to the vacuum, pair creation of electrons and positrons 
takes place and the insulation break down. The thresh- 
old for this phenomenon is known as Schwinger's critical 
field and is given by ii^th — "i^/e — 1.3 x lO^^V/cm. 
Since the critical field is extremely strong, direct exper- 
imental verification is still a challenge. On the other 
hand, in the condensed matter community, there is an in- 
creasing interest in non-equilibrium phase transitions and 
non-Unear transport in strongly correlated electron systems 
(Fig. [T]l Els] 131 El m 13. In the experiments, one also applies 
strong electric fields and the original insulating phase is de- 
stroyed. However, the threshold for dielectric breakdown 
is orders of magnitude smaller than the Schwinger mecha- 
nism in QED since the excitation gap is far smaller. This 
makes condensed matter systems to be an idealistic play- 
ground to test and develop theoretical ideas in non-linear 
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QFT. Non-linear physics has been studied rather indepen- 
dently in the two fields, high energy and condensed matter, 
during the past few decades, and several parallel ideas were 
developed. The aim of this article is to explain some of the 
correspondences (Table[T]l. 

PAIR CREATION IN STRONG ELECTRIC 
FIELDS 

Hesenberg-Euler' s effective action and the non- 
linear extension of the Berry's phase approach to 
polarization /El/-' 

We study lattice electrons in homogeneous electric 
fields. In the time-dependent gauge, this can be realized by 
adding a time dependent phase to the hopping term in the 
lattice Hamiltonian. For example, for a one-dimensional 
model, a typical Hamiltonian is given by 

i—l a 

i i 

We impose periodic boundary condition and the phase $ is 
proportional to the magnetic flux through the ring (i num- 
ber of sites). The time derivative of the magnetic flux is re- 
lated to the applied electric field through F(i) — eaE{t) = 
d^{t)/dt, where e is the charge quantum and a the lattice 
constant. U represents on-site Coulomb repulsion and Vi 
the local potential. The hopping term is set to unity. The 
Hubbard model ([/ > 0, = 0) at half-filling is in the 
Mott insulating phase for positive U in one dimension. 

Here, we study what happens to the an insulator when we 
apply strong electric fields. We denote the eigenstates of 
the Hamiltonian by |t/'„($)), ?i = 0, 1, . . . and study 
the time evolution starting from the groundstate | 'i/'o (*&))• 
The groundstate-to-groundstate amplitude defined by 

E{t) = (7^o(*(t))|e"''^° ''^*^'"''l^o(0))e'/o ^"(^(^"'^^ 

is of central importance. In the long time limit, an asymp- 
totic behavior (d is dimension) E{t) ^ e**^ is expected 
to take place where £ is the condensed matter version of 
the Heisenberg-Euler effective Lagrangian. The imagi- 
nary part describes quantum tunneling where T{F)/L'^ = 
2Im C{F) gives the speed of the exponential decay of the 
vacuum (groundstate). This quantity is proportional to the 
decay rate of the Loschmidt Echo L{t) — The 
real part Re£ is written in terms of a non-adiabatic phase 



Table 1 : Related ideas in strong field physics 



High Energy 

Schwinger mechanism in QED 
Heisenberg-Euler effective Lagrangian 
Vacuum polarization 
Pair creation in interacting systems (e.g. QCD) 

Dirac particles in circularly polarized light 
Furry picture 



Condensed Matter 

Landau-Zener tunneling in band insulators 
Non-adiabatic geometric phase, Loschmidt Echo 

Extended Berry's phase theory of polarization 
Many-body Schwinger-Landau-Zener mechanism 
in strongly correlated system 
Photovoltaic Hall effect 
Floquet picture 



called the Aharonov-Anandan phase (which we denote 7) 
that the wave function acquires during the time-evolution. 
For band insulators (U — 0) in dc-electric fields, the effec- 
tive Lagrangian becomes Q 

f dk j(k) 
Re£(F = -F 7^^, (3) 
Jbz {^T^r 27r 
r fiu 1 

Im£(F) ^ -^^y^^^-ln[l-p(fe)], (4) 

where the momentum integral is over the Brillouin Zone 
(BZ). There is an interesting parallel theory developed in 
the condensed matter community. This is known as the 
Berry's phase theory of polarization|[8l|9l[l0l[TTl[l2l, where 
the ground-state expectation value of the twist operator 
e~*T7^^ which shifts the phase of electron wave functions 
on site j by — ^j, plays a crucial role. It was revealed that 
the real part of a quantity 

^^_!ln(0|e'^T*|0) (5) 

ZTT 

gives the electric polarization Pel = — Rew IJOJ while itfs 
imaginary part gives a criterion for metal-insulator transi- 
tion, i.e., D = 47rlmw is finite in insulators and divergent 
in metals 111 11 . The present effective Lagrangian can be 
regarded as a non-adiabatic (finite electric field) extension 
of w. To give a more accurate argument, we rewrite the 
effective Lagrangian in the time-independent gauge 

C{F) ^ -^ln((0|e-^-(^+^^)|0)e^-^'') (6) 

for d = 1. Let us set r = h/LF and consider the 
small F limit. For insulators we can replace H with the 
groundstate energy Eq to have C{F) ^ wF in the linear- 
response regime. Thus the real part of Heisenberg-Euler's 
expressionflJl for the non-linear polarization P{F) — 
—dC{F)/dF naturally reduces to the Berry's phase for- 
mula Pel in the small field limit F ^ Q. Its imaginary part 
gives the criterion for photo-induced metal-insulator tran- 
sition, originally proposed for the zero field case. 

Many-body Schwinger-Landau-Zener mecha- 
nism in strongly correlated insulators: 

Next, let us consider dielectric breakdown in a strongly 
correlated system. In the one-dimensional Mott insulator 



where the groundstate is a state with one electron per site, 
the relevant charge excitations are doublons, i.e,. doubly 
occupied sites, and holes, i.e., sites with no electron. Pairs 
of doublons and holes play a similar role as the pair of 
electrons and positrons in the Schwinger mechanism. In- 
deed, it has been shown that dielectric breakdown in Mott 
insulators takes place due to pair production of charge ex- 
citations through quantum tunneling, which is called the 




Figure 2: (a) Many-body energy levels against the com- 
plex AB flux <i> for a finite, half-filled ID Hubbard model 
(L = 10, N-^ = = b, U = 0.5). Only charge exci- 
tations are plotted. Quantum tunneling occurs between the 
groundstate (labeled as n = 0) and a low-lying excited state 
(n = 1) as the flux — Ft increases on the real axis, 
while the tunneling is absent for the states plotted as dashed 
lines. The wavy lines starting from the singular points (x) 
at ^{t*) represent the branch cuts for different Riemann 
surfaces, along which the solutions n = and n — 1 
are connected. In the DDP approach, the tunneling fac- 
tor is calculated from the dynamical phase associated with 
adiabatic time evolution (DDP path) that encircles a gap- 
closing point at ^{t*) on the complex $ plane, (b) Thresh- 
old electric field obtained by the imaginary time method 
(solid) and the naive Landau-Zener formula (dashed). 



many-body Schwinger-Landau-Zener mechanism f3','4','5l. 
The tunnehng probabiHty for the lowest excited state in the 
one-dimensional Hubbard model was recently calculated 
by combining the imaginary time method with the exact 
Bethe ansatz solution |5|. 

The imaginary time method has been widely used in 
quantum physics, e.g., in quantum chemistry lfT4l [Tsl and 
in high energy llT6l[T7iri8] . In high energy, the tunneling 
rate for Dirac particles in ac-electric fields was calculated 
by analytically extending the time-evolution to complex 
time processes. This method can be readily applied to the 
one-dimensional Hubbard model since the wave functions, 
including the excited states, are analytically known. Using 
Bethe ansatz, these states can be described as a linear com- 
bination of plane waves with the momentum (rapidity) de- 
termined by the Lieb-Wu equation 1 19 1. In the infinite size 
hmit, the Lieb-Wu equation can be solved even in the pres- 
ence of a complex phase, i.e., $ ^ $ + i^* 1201 |5l. In the 
half-filled Hubbard model, the lowest energy excited state 
is given by the string solution ET\ l22l 123. ,24j . We denote 
the energy difference between this state and the ground- 
state by A = El — Ed- For dc-electric fields, quantum 
tunneling to the lowest energy state dominates and one can 
express the decay rate by 

r/i ^ _ ^ ln[l - exp(-^^^th/i^)] , (7) 

where Fth is the threshold field for this process and a a nu- 
merical factor. In the imaginary time method, an adiabatic 
evolution to a gap closing point in the complex ^ + W 
space is performed and the tunneling probability depends 
on the dynamical phase. The gap closing in complex $ has 
been studied in ref. f20l, while in ref. f25l a related problem 
was studied in order to calculate the localization length in a 
Mott insulator. The Bethe ansatz solution for the Hubbard 
model with finite ^ is represented by a rapidity distribution 
with end points ±a + ib. The imaginary part b is related to 
^ hy a condition that the charge quantum number should 
remain real. The threshold for the tunneling is given bytSJ 
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where A is the charge gap (Mott gap) f2E\, and v is the 
slope of the adiabatic levels {v ^ 2 when U is small and the 
system size is small). From Fig.|2](b), one notice that the 
result of the imaginary time method and the Landau-Zener 
formula coincides in the small U limit. However, when the 
interaction becomes stronger, large deviation takes place, 
with different asymptotic behaviors F^"^^ oc U — const, 
and F^ cx [/^. Dielectric breakdown in the Hubbard 
model has been studied numerically and threshold behav- 
iors in the decay rate |4| as well as in the current 1271 161 are 
observed. Further research to explore the strong U limit is 
yet to be done. 

In ac-electric fields, the tunneling behavior can still be 
realized. The Keldysh line in Fig.[T]corresponds to the line 
where the adiabaticity parameter |>28| 7 = ^/F^ crosses 
unity (^ is the localization length in the present case). For 
7^1, the tunneling process is the dominant cause of car- 
rier production, while when 7^1, multiphoton absorp- 
tion wins. In the next section, another interesting effect of 
the ac-electric field is discussed where light can be used to 
control the topology of the band structure. 

PHOTOVOLTAIC HALL EFFECT 
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Its tZ-dependence is plotted e in Fig. [2] (b) (solid line). 
This indicates a collective nature of the breakdown (i.e., 
the threshold much smaller than a naive U). In other words, 
the tunneling takes place not between neighboring sites, but 
over an extended region due to a leakage of the many -body 
wave function, where the size is roughly given by the local- 
ization length [25|. One can compare this result with the 
naive estimate for the tunneling threshold obtained by the 
Landau-Zener formula ^ I?) 

(A/2)2 



(9) 
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Figure 3: Fig. 5 (upper) Greenfs function in the Floquet 
picture, (lower left) Diagram for conductivity, (lower 
right) Photo-induced Berry curvature (Chern density) of 
graphene in circularly polarized light. 

Dirac electrons are now becoming one of the central top- 
ics in condensed matter after its experimental realization 
in graphene ||29l . Electrons in strong background electric 
fields is an important problem and can be experimentally 
studied by optical and transport methods. In non-linear 
QED, this problem has been studied by Volkov for ac- 
fields where he found that electron wave functions acquire 
a non-trivial phase[30]. However, in condensed matter, the 
electrons subject to strong lasers behave differently from 
the Volkov state of relativistic Dirac particles. This is be- 
cause the velocity, or the "speed of light", of Dirac elec- 
trons in condensed matter systems are much slower com- 
pared to the actual speed of light and the ac-field breaks 
the "Lorentz symmetry" of the electron. This leads to a dy- 
namical gap opening in the pseudo-energy diagram lISTl lTl. 
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Figure 4: (upper) DC current. (lower left) IV- 
characteristics. (lower right) Conductance. 

Especially, in a circularly polarized light, a gap opens at 
the Dirac point ||7|- This has an important physical conse- 
quence since a gap of a 2h-1 dimensional Dirac electron is 
related to parity anomaly and is detectable through trans- 
port measurements, i.e., the Hall effect. In 2+1 dimension, 
the Hall conductivity can be written as a momentum in- 
tegral of the Berry curvature (~ Chern density) over the 
Brillouin zone. This is known as the TKNN formula|32|, 
and is know extended to ac-driven transport via the Floquet 
picture (^ Furry picture) |7| 
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Uk) [Vfc x^„(fc)]^. (10) 



Here, Aa{k) = -i(($Q(fc)|V^|<i>Q(fc))) is the photo- 
induced artificial gauge field. In the Floquet picture, the 
Green's function incorporates the effect of photon absorp- 
tion and emission (Fig. [3] (a)), and Hall conductivity is 
given by the bubble diagram in the non-interacting case, 
which is nothing but the parity anomaly diagram. The 
photo-induced Berry curvature shown in Fig.[3](c) acts as 
an artificial magnetic field and becomes finite when the cir- 
cularly polarized light is introduced. 

The current in the presence of circularly polarized light 
in a graphene ribbon attached to two electrodes is plotted 
in Fig. [4] The calculation has been done by combining the 
Keldysh green's function method with the Floquet picture. 
The Hall current, which is originally absent, increases as 
the strength of light becomes stronger. The numerical re- 
sult supports our understanding of the photovoltaic Hall ef- 



fect obtain by the extended TKNN formula (eqn.( 10 1) 
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